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ABSTRACT 

The acceleration of the expansion of the universe, ascribed to a dark energy, is one of the most 
intriguing discoveries in science. In addition to precise, systematics controlled data, clear, robust 
interpretation of the observations is required to reveal the nature of dark energy. Even for the simplest 
question: is the data consistent with the cosmological constant? there are important subtleties in the 
reconstruction of the dark energy properties. We discuss the roles of analysis both in terms of the 
Hubble expansion rate or dark energy density pde{z) and in terms of the dark energy equation of 
state w{z)^ arguing that each has its carefully defined place. Fitting the density is best for learning 
about the density, but using it to probe the equation of state can lead to instability and bias. 
Subject headings: dark energy — cosmologyiobservations — cosmology:theory 



1. INTRODUCTION 

The acceleration of the expansion of the universe rep- 
resents a challenge to our understanding of fundamental 
physics. Whether the resolution of this mystery lies in 
high energy physics, the theory of gravitation, the nature 
of the quantum vacuum, extra dimensions, etc. remains 
unknown. The simplest m odel, and possibly one mo- 
tivated from string theory l)Kachru 2 003). is Einstein's 
cosmological constant, though it too is fraught with com- 
plications, e.g. the fine tuning and coincidence problems 
ilWeinberg lOMlCarroll 20nl . 

Nevertheless, the cosmological constant serves as a 
benchmark, a first indication for how drastically we 
might need to extend our physics framework to under- 
stand the acceleration of the universe. Additionally, be- 
cause this model offers such definite and unvarying pre- 
dictions for the dark energy properties, it gives a robust 
target. In this paper we consider various aspects of in- 
terpretation of the data, present and future, that aid or 
obscure the goal of determining the nature of dark en- 
ergy and its consistency with the cosmological constant. 
In the sense that conflicting or biased interpretations of 
the data, rather than the data itself, affect what we de- 
termine to be the reality, we refer to this problem as 
"deconstruction" of dark energy. 

2. RECONSTRUCTION AND DECONSTRUCTION 

The cosmological constant has an energy density that 
is constant in time, p{z) = p(0) in terms of redshift z, a 
pressure the negative of the energy density, p = — p, or 
equivalently an equation of state ratio (EOS) w = p/ p = 
— 1 constant in time, and no spatial fluctuations in p or 
w. So we could think of testing whether the dark energy 
is a cosmological constant by seeing if the observations 
are consistent with either a constant density or w = —1. 



But this obscures the physics behind the dark energy. 
Even if we are so modest and undemanding as to only 
ask if Wconst = — 1 or not, we can be fooled by time 
varying models mimicking this value. 

We illustrate this in Figure ^ The two dark en- 
ergy models shown, one a high energy physics model 
parametrized with a time evolution w{a) = wq -I- Wa(l — 
a), where a = 1/(1 + z) is the scale factor of the universe, 
with (wo, Wa) — (1.5, —1.5), a nd one an extended gravity 
model (case 3 of lLinder 2004() modifying the Friedmann 
expansion equation, both possess a strong time varia- 
tion in the equation of state. However, in comparison to 
present data both appear consistent with a cosmological 
constant when interpreted in terms of a Wconst • 
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2.1. Constant equation of state Wconst = —1 

The first point to make is that fitting the data in terms 
of a constant EOS parameter is insufficient. Many mod- 
els with time variation w' can be fit, for data of limited 
precision or over a limited redshift range, by a Wconst- 
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Fig. 1. — Dark energy models possessing strong time varia- 
tion can still appear consistent with cosmological constant when 
interpreted only in terms of a constant equation of state. The 
two models shown here are consistent with all current data and a 
'"'const = — 1 cosmological constant model. 
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Both models are consistent, to 0.1%, with the distance 
to the last scattering surface of the cosmic microwave 
background (CMB) in the cosmological constant model 
with the same matter density. Even the n ext genera- 
tion Planck experiment ijsee Planck websita l will not be 
able to distinguish these models from the cosmological 
constant on this basis. With regard to large scale struc- 
ture, the linear theory growth factors agree to better than 
4%. This roughly corresponds to requiring knowledge of 
the mass fluctuation amplitude erg a-t this level, but the 
current systematic uncertainties appear to be of order 
20%. Most constraining are the Type la supernova dis- 
tances, for which the time varying models deviate from 
the cosmological constant by under 0.1 magnitude. This 
is currently at the la limit of precision. 

So if our universe followed either of these dark energy 
models (or a host of others that could be chosen) then 
by using the Wconst paradigm we would happily interpret 
the data as pointing to a cosmological constant and miss 
the key clues to new physics. 

Of course as cosmological probes become more pre- 
cise, we will be able to distinguish the illustrated models 
from the cosmological constant. However if the interpre- 
tation does not include the possibility of time variation 
in the dark energy equation of state we will again miss 
the physics and think the particular models shown are 
Wconst = —1.2 models. And indeed dark energy mod- 
els with a less extreme time variation than shown in 
Fig. n might well look like Wconst = — 1 to within the 
data precision. Moreover, many models have an asymp- 
totic or attractor behavior that brings the EOS very 
near the cosmological constant value for a period of the 
expansion history. Exampl es includ e the li near poten- 
tial mo del llLinde 1987t iKal losh et al. 200,"^ . the cyclic 
model (jSteinhardt fc Turok 2004^) . the "ripstop" model 
where particle production offset s the superac celeration 
of phantom energy with w < — 1 (|Li nder 20 (M). a nd cer- 
tain scalar-tensor theories (jMatarrese et al. 2004J) . 

For a window onto the new physics we need to look 
for both the present or averaged EOS value and a mea- 
sure of its time variation, say w' . Note that throughout 
this paper we refer to the EOS function w{z) as an ef- 
fective EOS, defined in terms of the expansion history 
a{t), and not necessarily rest ricted to a physical coin po- 
nent lik e a scalar field. From lLmder fc Jenkins 20031 (cf. 
ISaini et al . 2003) we have the effective EOS defined by 



Note that this is closely related to the inflationary theory 
measure of time variation: 



1 dhiSH'^ 
w{z) = -1 + -- 



(1) 



3<iln(l + z)' 
where the Friedmann expansion equation reads 

H^{z)/H^^nM{l + zf + SH\z). (2) 

The Hubble parameter H = d/a and Qm is the matter 
density, so SH^ represents our ignorance of either compo- 
nents beyond matter or of extensions to general relativity 
that modify the Friedmann equation. 

A very robust and convenient measure of the time vari- 
ation of th e dark energy EOS is given by the two param- 
eter form l|Linder 2003tlLinder 200^ 

w{z) = Wq + Wa{l - a) = Wo + WaZ/{l + z) . (3) 

We see that this implies a characteristic variation 

Onw = -"71 = awa- (4) 

a In a 



Onw — ~dw/dN, 



(5) 



where N = din a counts the number of e- folds of ex- 
pansion during inflation. Just as in inflation, there is 
"running" ~ a variation of parameters during the accel- 
eration - and we need to choose when to evaluate the 
derivatives, e.g. at A^ = 60. The derivative of the EOS 
will also be a function of the redshift when it is evaluated, 
and for simplicity we define a number from the function: 
w' = Onw^z = 1) = Wa/2. This is motivated by z « 1 
being approximately when dark energy begins to become 
dynamically important for the expansion. 

So we have argued that to learn the direction of new 
physics we need to consider not only the value of Wconst 
but not deny the possibility of time variation w' . That is, 
if we want to say that observations point to the cosmolog- 
ical constant, we need to see that w' ~ 0. A much more 
modest goal is merely to ask whether the observations are 
consistent with the cosmological constant. In this case, 
we might consider looking for whether the dark energy 
density is constant, rather than addressing the physics 
question of its equation of state. 

2.2. Constant density 

Note that reconstructing a constant energy density is 
only a necessary, not sufficient, condition for a cosmologi- 
cal constant. Some models might arrange an effective en- 
ergy density nearly constant with time yet with effective 
pressure not equal to the negative of the energy density; 
this holds in some k-essence and tachyon field models (cf. 
iScherrer 20041) . So mapping the density alone contains 
insufficient information if we want to remain open to the 
possibility that the unknown physics could involve com- 
ponents other than canonical, minimally coupled scalar 
fields or involve extensions to gravity. 

But in asking for consistency, not insight, seeking 
to rec onstruct the dark ener gy density behavior (see 
e.Ef. iWanEf fc Tesfmark 2004t iWang fc Freese 20041: 
iHuterer fc Coorav 2004(1 is a perfectly acceptable 
method. Testing for the EOS w(z) also allows one to 
check consistency (plus more of course), and one can 
also build up the density history through an integration: 



p{z) = piO)e'Jo 



3fg'dz'[l+wiz')]/{l+z') 



(6) 



Of course this will diminish the sensitivity of the re- 
construction of p(z) som ewhat, though not drastically 
(|Huterer fc Coorav 2004|) . Also, it may be easier to test 
for a specific v alue, e.g. w = — 1, than for a property, e.g. 
p = constant l)Huterer 2004(1 . Generally, if one is only 
interested in whether the density is constant, use recon- 
struction in terms of the density (though in terms of EOS 
is not unreasonable). However, if one is interested in the 
physics behind the dark energy, the EOS reconstruction 
is strongly preferred. 

This is because going from p{z) to w{z) requires dif- 
ferentiation, whether explicitly or implicitly, and this 
leads to instability and bias in the results. This is 
what we referred to in the introduction as the decon- 
struction problem. The dange rs of explicit differentia- 
tion of the data are well k nown J Huterer fc Turner 1999t 
iWeller fc Albrecht 20021: iTeemark 2002() . but those of 
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implicit differentiation hav e been illustrated more re- 
cently. I.Tonsson et al. 200l pointed out that assuming 
a functional form for the density, or Hubble parameter, 
could distort conclusions regarding the equation of state, 
even apparently preferring a time evolving dark energy 
to fit data of a cosmological constant universe. 

Here we extend this to the general translation from a 
density reconstruction to an equation of state behavior. 
Consider a set of values {pi}, representing the density in 
redshift bins centered at z^. To obtain the EOS, one can 
imagine fitting the reconstruct ed pj with a polynomia l 
or spline interpolation p{z) (cl. IWang fc Tegmark 200^) 
and then carrying out the derivative to form 



1 / sdp 



In terms of the binned values, 



w{z) = -1 + 



1 + Z J- Pi+1 

3p.i Iz^+i 



(7) 



(8) 



Zi+l 



-Pi 



3(zi+i - zf 






(Zj+l - ZiY 

■i{z-Zif 

(Zi+1 - ZiY 
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where the second derivatives are determined in terms of 
the free parameters (bin values pi) through the cubic 
spline formula and boundary conditions. Thus we can 
write Eq. |H1 illustratively as 



w(z) 



-1 



[l/(3p)](l + z)5]c, 



Pi 



(9) 



where pi are the free parameters of the problem and ci are 
weight functions. One can imagine these being more gen- 
eral than cubic spline coefhcients and instead being either 
polynomial expansions in redshift, e.g. q = hi{z — Zi)"-, 
or involving orthogonal basis functions, so the following 
conclusions should be fairly general. ^^_^^_^ 

Following the method of Jonsson e t al. 20(M now con- 
sider the stability of this deconstruction to small pertur- 
bations 6 Pi . The quantity in braces in Eq. |S1 (or its gen- 
eralizations) will involve terms like Spi, {dp'l / dpi)5pi ^ 
6 Pi. However, formally, because of the 1 + z factor in 
front of the braces, a small perturbation to the fit to 
the data in terms of pi will run away at high redshift in 
terms of w(z), showing the instability of this approach^. 
(Note that redshifts not mu ch higher than 1 may be suf- 
ficient; iJonsson et al. 20M find instability regions where 
\w\ > 15 for z < 1.8.) 

Furthermore, this approach not only mediates against 
the cosmological constant (or any Wconst model), but fa- 
vors a model where the EOS crosses from less than —1 to 
greater than — 1 : precisely what is seen in current recon- 
structions. This occurs because the quantity in braces is 
redshift dependent and will have a zero at some redshift 
z,. That is, for generic sets {pi} one will find w(z*) = —1 
(not necessarily within the redshift range of the data). 
Thus the deconstruction legislates for crossing the value 

^ It also may be problematic to take a linear interpolation as 
one of the spline boundary conditions as Wang & Teemark 2004 
do. (Note that the p'{z = 0) = p{zi)/zi appearing there is a 
typo for p'{z = 0) = [p(2i) — p(0)]/zi.) This may allow statistical 
excursions of the density off the cosmological constant behavior to 
cause w yt —1. 



— 1. Combining these two effects together, a small per- 
turbation 6pi will generate a larger 6w, forcing w off from 

— 1 (if the cosmological constant were the true fit) and 
forcing w to cross —1 - both leading to evolution. 

In summary, using pi fitting to find p{z) is fine, but 
using Pi to discuss w(z) is perilous. The latter generates 
errors of two types: false negatives where a true cosmo- 
logical constant appears not to be so, and false positives 
where any model looks, over some finite redshift interval, 
not very distinct from a cosmological constant! 

That said, there appears to be a simple, if partial fix. 
We can rewrite Eq.Qin terms of a logarithmic derivative 
with respect to scale factor a, rather than 1 + z. If bins 
are defined in terms of p{a), then the high redshift secu- 
lar instability does not appear, and this should mitigate 
any tendency to cross the value —1. General issues of 
numerical differentiation, bin to bin, instability and the 
need for careful treatment of spline boundary conditions 
remain, but are possibly less severe. 

2.3. Density or Equation of State? 

What about the breadth of models encompassed by 
the tw o approaches? It has sometimes been claimed 
IjWang & Teemark 2004) that w{z) suffers from not cov- 
ering more esoteric possibilities like negative density. 
This is not wholly true, and even so is more of a fea- 
ture than a bug. From Eq. ^ one sees that negative 
densities can be handled straightforwardly in terms of 
w(z). For example, consider a model where we (mistak- 
enly) think riM = 1.2 and find the best fit to data gives 
SH^ — —0.2(1 -t- z)"^. Clearly this looks like a negative 
density. But [(1 + z) / S H^]dS H^ / d{l + z) is well defined 
and gives w{z) = 0: exactly as expected for this extra 
(negative) matter component. 

Where the w{z) parametrization does blow up is not 
for negative density but for density crossing thro ugh zero. 
In the parametrization of iWang fc Tegmark 2 004. they 
normalize the density by the present value, using X{z) = 
p{z)/p{0), with of course X(0) = 1. In the negative 
density illustration above, X will then always be positive. 
It is only when X goes from its defined value of 1 at z = 
to at some redshift z, that wlz^) blows up to infinity. 
But we can argue this is a feature not a bug. 

For example, in th e linear potential model analyzed 
in lKallosh et al. 200,'^ the dark energy density hits zero 
very shortly in cosmic time before the final "cosmic 
doomsday" collapse. By parametrizing the w{z) given 
by the linear potential (e.g. in terms of wo and Wa, con- 
trary to the claim by IWang fc Tegmark 2 004 that this 
model cannot treat the collapse case) and asking when 
it blows up (in the future), one obtains a very accurate 
estimate for the doomsday time. 

Furthermore, having p cross through zero is not merely 
mildly esoteric, but violates either the Big Bang condi- 
tion or the continuity equation 

dp/d\n{l + z) = 3{p + p). (10) 

Dividing through by p(0), one has (iX/dln(H-z) = 3{X+ 
p/p{0)). For acceleration, we take the pressure p to be 
negative and so at the time when X = 0, the left hand 
side will be negative and p will be driven further negative 
at high redshift. This negative energy density will cause 
the expansion rate H to go to zero at some point in 
the past, describing a bounce universe rather than a Big 
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Bang model. If one insists on allowing the density to 
cross through zero then the price is either abandoning 
the continuity equation or the Big Bang early universe. 
And if one abandons the continuity equation then the 
density alone is insufficient since pressure is unknown, 
leading back to an equation of state formulation. 

Finally, one can argue that the Friedmann equation of 
motion explicitly involves both energy density p and pres- 
sure p, so a parametrization w{z) = p{z)/p{z) is closer 
to the physics. An example of this was pointed out at 
the beginning of ii2.2l where a constant energy density 
may not necessarily be matched by a constant negative 
pressure of equal magnitude. Fundamentally, the passive 
gravitational energy p+p and the active gravitational en- 
ergy p + 3p both explicitly depend on more than merely 
the density, they involve an equation of state. 

3. CONCLUSION 

Understanding, or at least obtaining insight into, the 
nature of dark energy will be the great challenge of 
physics in the next decade. We must be sure that we 
ask questions in such a way that the answers we derive 
are not deconstructions - subjective interpretations - but 
faithful reconstructions of aspects of the true physics. 
This includes allowing explicitly for the possibility of 
time variation in the equation of state of the dark en- 
ergy. Using a Wconst will only tell us whether the data 
are consistent with the cosmological constant, not teach 
us that another, mimicking model is not the true answer. 



If our main concern is whether the dark energy 
density is constant in time, we can use reconstruc- 
tion of the density values in (preferably uncorrelated) 
bins in expansion factor or a principal component 
analysis of the density (see iHuterer fc Coorav 20041: 
iHuterer fc Starkman 200.1() . Using the equation of state 
w{z) in either bins or eigenmodes will give almost as ac- 
curate results, plus additional information. Crosscheck- 
ing one versus the other may also test exotic models in- 
volving violation of the continuity equation. 

In our quest to understand new physics, key clues are 
carried by the equation of state. To reveal this one 
should fit for w{z) directly. Using other, intermediary 
parametrizations such as density on the way to the equa- 
tion of state can lead to unstable and biased solutions, 
both type 1 (false positive) and type 2 (false negative) er- 
rors with respect to determining if the data are consistent 
with no time evolution, w{z) = — 1, like for the cosmo- 
logical constant model. Of course systematics errors in 
the data themselves can also lead to false evolution. 
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